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We develop a systematic approach to calculating the electrostatic force between point charges 
in an arbitrary geometry with arbitrary boundary conditions. When the boundary is present, 
the simple expression for the force acting on a charge as "the charge times the field it is placed 
in" becomes ill-defined. However, this rule can be salvaged if the field in question is redefined to 
include all the terms that do not diverge at the charge position, in particular, those due to the 
charge itself. The proof requires handling the self-action energy divergence for point charges, 
which is accomplished by means of a geometrical regularization. 
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I. INTRODUCTION 

It is trivial to determine the force exerted by an exter- 
nal field 1 on a point charge in an otherwise empty space: 
by definition, "the force is equal to the charge times the 
field it is placed in" . In particular, if the field in question 
is created by some other point charges, this rule, known 
by many from high school, still holds. 

However, the situation changes drastically when a set 
of point charges creates the field inside an arbitrary do- 
main with a boundary of some physical origin (reflected 
in the appropriate boundary conditions). Now the very 
notion of "the field the charge is placed in" becomes not 
well defined. For example, a naive treatment of a single 
charge problem might lead one to an entirely wrong con- 
clusion that, since all the field in the problem is due to 
the charge itself (there are no other sources!), "the field 
it is placed in" is zero, so there is no force at all. 

A slightly more sophisticated physicist would argue 
that only the part of the field which diverges as 1/r 2 
near the charge is really created by it, while the rest is 
due to the boundary conditions, which represent math- 
ematically the rearrangement of other physical charges 
at the boundary. Therefore it is precisely what remains 
after subtracting the singular part that now gives "the 
field the charge is placed in" . Unfortunately, such treat- 
ment leaves one in a somewhat awkward position of, first, 
calculating potentials and fields rigorously, and then low- 



ering the plank and using hand-waving arguments to de- 
rive forces from them. It is also not clear whether the 
guess about which part of the total field contributes to 
the force is always valid. 

Thus it seems appealing to show that the physical 
arguments can be backed by an accurate mathematical 
proof demonstrating that the adjusted rule, "the force is 
equal to the charge times the part of the field that does 
not diverge at the charge's location" , is either universal 
or limited by certain conditions. To do this, one should 
turn to the most fundamental energy conservation argu- 
ment which gives the force as the negative gradient of 
the energy in the charge's position. This approach also 
does not turn out to be straightforward, since the energy 
is infinite in the presence of point charges due to their 
self-action. 

Perhaps because of these difficulties, as well as of a 
misleading apparent simplicity of the problem, our liter- 
ature search, which encompassed, in particular, Refs. ^, 
|, [| |, |§ §, ]| g g, [fi] and many other books on the sub- 
ject, revealed no ready result (except for one small piece 
in Ref. ^] which we discuss in Sec. IV). So we give a care- 
ful derivation of the general expression for the force on 
point charges in this paper. It consists of a regularization 
of the problem, calculation of the force from the (regu- 
larized and finite) energy, and then taking the singular 
limit. The result agrees with one's intuitive expectations. 
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Hn fact, the word "external" in any case means that the sources 

of the field, such as some boundaries, are far away from the charge. 



II. ELECTROSTATICS PROBLEM WITH 
VOLUME POINT CHARGES: POTENTIAL 
AND ENERGY 

Consider an arbitrary 3-dimensional domain D with 
the perfectly conducting boundary S and some N point 
electrical charges inside. The electrical potential, "0(r), in 
this case is determined by the following Dirichlet bound- 
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ary value problem (we use SI units throughout the pa- 
per): 



1 N 

Aip = } qi 6 (r - n) , r, r, e 



D: 



4> 



= 0. 



(1) 

(2) 



Here r = xe x + ye y + ze z is the vector radius of a point, 
and = X{G X + yie y + z;e 2 specifies the ith charge po- 
sition, with e a , a = x,y, z, being the unit vectors in the 
direction of the corresponding Cartesian axes. 

By superposition principle, the potential f/ , ( r ) is j us t 
the sum of the potentials induced by each charge sepa- 
rately, 
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where /c = l/47T£o, and G R (r, Tj) is the regular part of the 
Green's function G(r,Tj) of the corresponding boundary 
value problem [set Qj = 1, = 0, i ^ j in Eq. ([!])]. Both 
functions are, of course, symmetric in their arguments, 

G(r,r J ) = G(r J ,r), G R (r,rj) = G*(r,-,r). (5) 

Furthermore, we can rewrite Eq. (|]) splitting the poten- 
tial in a sum of its singular and regular parts, 
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where iPr(v) is a regular function satisfying the Laplace 
equation everywhere in D [by continuity, this holds also 
at any regular point 2 of the boundary S, although this is 
irrelevant to our discussion]. Note that both the poten- 
tial, ip, and its regular part, ipR, depend actually on the 
positions of the charges as well as on the observation 
point r, which is reflected in the full notation, 
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We assume that the potential is known and are inter- 
ested in finding the force F* acting on the charge qi. From 



the energy conservation for the considered problem, the 
force is given by (cf. Ref. ||) 3 
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F l = -^-W D , tt- = ^~e x + — e y + —e z , (10) 
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where Wr> is the energy of the field in the volume D, 



* J D 
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Note that we alternatively write V or d/dr for the gra- 
dient, whatever seems proper in a particular expression. 

The problem is, however, that the above integral obvi- 
ously diverges due to self-interaction of the point charges 
(the energy of a single point charge is infinite). We are 
going to show that even though the energy for a given 
point charge distribution is infinite, the difference be- 
tween its two values corresponding to two different but 
close charge configurations is finite for any charge con- 
figuration and boundary shape, and the force is also finite 
due to that, in accordance with common intuition. 



III. REGULARIZED ENERGY AND THE 
FORCE ON THE CHARGES 

We surround each volume charge qi by a small sphere 
Sf of radius e; we write Df for the ball inside it. We define 
D e as D without all domains D\ , and S e as a union of S 
and all spherical surfaces S% (see Fig. [j]). In effect, S e is 
the boundary of D e . 

Using Eq. (10), we may now define the force acting on 
the charge qi as 



F< = lim Ft = -Um-J?-Wf» 



(12) 



where is the regularized energy, that is, the energy of 
the field in D c , which is finite. It is important to note the 




FIG. 1: Volumes, surfaces, and normal directions involved. 



2 We allow for the boundary singularities, such as sharp edges 
and sp ikes provided that the Meixner type finite energy condi- 
tion [hl| is satisfied near them; in particular, the domain D can be 
infinite. 



3 ft is important to understand that Eq. (|l(j) is a definition of 
a mathematical object that we would like to correspond to the 
physical force. If there are no external fields [Eq. (0) ensures this 
in ourdiscussion] , it will turn out that the conclusions derived from 
Eq. (hoh are physically meaningful and validate the definition. 
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order of operations in Eq. ([L2|) : first take the gradient of 
the regularized energy in the charge position, then take 
the (singular) limit. In principle, we also have to show 
that the final result does not depend on the regularization 
chosen, but this task is not easy. We will return to it 
briefly later in this paper. 

In view of Eq. (pT ) and the fact that the total po- 
tential given by Eq. (|) or Eq. (||) is regular in D e , the 
regularized energy is 
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ipAipdV; (14) 



n is the direction of the outward normal to S e (and thus 
the inward normal to the spheres For an infinite 

domain D it is assumed here that the potential and its 
gradient drop at infinity fast enough to make the con- 
tribution of integrating over the sphere of a large radius 
vanishing in the limit, which assumption has to be veri- 



fied in each particular case. 

Since ip is harmonic everywhere in £) e , the volume in- 
tegral on the right of the previous equality vanishes; the 
remaining surface one is represented as 
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and then, because of the boundary condition, Eq. 
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as 
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We are ultimately interested in the limit e — > 0, so we 
need to calculate only the quantities which do not vanish 
in this limit. The area of integration in each term of the 
above sum is 0(e 2 ), therefore we need to keep track of 
the integrands that grow at least quadratically in e _1 . 
Bearing this in mind and using Eq. (^|) for the potential, 
we can write the surface integral in Eq. ( |l6| ) as 
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The first term in the above expression is, in fact, a reg- 
ularized self-energy of the fe-th charge, sel f. Doing 
an elementary integration, we immediately find that 
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The only feature of the regularized self-energy given by 
Eq. dig) important for our derivation is that it does not 



depend on the position of the charge q k , i. e., on the 
vector radius r k . 

The second term of the r. h. s. of Eq. ( |l7| ) can also 
be simplified if one notices that both ipn and 1/ 1 r — rj\, 
j ^ k, are regular on S^. and in D e k . Therefore their 
change within the small surface Sf. is of order e. Thus 
Eq. (|l7|) may be rewritten as 
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and the integration here yielding the factor An is again an elementary one. This asymptotic equality may be differ- 
entiated in Vi with the same estimate of the remainder term. 

Introducing now the last expression into the Eq. (jig), we obtain: 
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Equation (|20|), in its turn, is inserted in the Eq. (|12|) for the force; as shown, the self-energies do not depend on 
the charge positions, hence, although diverging in the limit e — > 0, they do not contribute to the force. The rest 
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is pretty straightforward, except one has to be careful when differentiating the last term on the right of Eq. ( pQ ) 
with k — i: as seen from Eq. (g|), in this case r 4 stands for two (and not one!) arguments of ipR, namely, Vfi( r i) = 
^>ii(rj, ri, . . . , Tj, . . . ,rjv), and both of them have to be differentiated. Bearing this in mind, the expression for the 
force finally becomes: 
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This is the general result for the electrostatics which can be transformed further in some nice way. Indeed, the 
direct substitution of the expression for ip R from the Eq. (||) in the the Eq. (|l]) provides the force in the form 
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and we have used the symmetry property of Eq. (^) to obtain the second equality here. To make the result even more 
physically transparent, we rewrite Eq. (E2), in its turn, in the following way: 
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Note that the last expression, indeed, coincides with our 
intuitive guess for the form of the force. 



IV. DISCUSSION 

Our first remark on the expressions for the force in 
Eqs. (21)~(|23]) is that for the charges in a free space (vol- 



ume D is the whole space, no boundaries present) appar- 
ently Gi?(r,rfc) = 0, i)jr = 0, and the classical Coulomb 
formula for the force is restored. 

Next, Eq. ( p3| ) shows that the rule "the force is the 
charge times the field it is placed in" does work if one 
counts the regular part of the field produced by the charge 
in question as a part of the "field the charge is placed in". 
It also makes up to some "minimal principle" , namely: 
to get the right answer for the force, one should throw 
out of the field only the part which otherwise makes the 
result infinite, and nothing beyond that. As we mentioned 
in the Introduction, this result is supported by physical 
intuition. It becomes even more transparent if one notes 
that the singular part of the field thrown out is radial, 
and the radial field produces no force. 

A contribution of the regular part of the field created 
by a charge to the force acting on it is especially impor- 
tant in the case of a single charge, as one may see from 
the simplest example of a charge near a conducting plane. 



It is exactly the regular part of the field produced by the 
charge in question (equal to the field of the image charge) 
that gives the whole answer when no other charges are 
present. 

Finally, an important question is how robust our regu- 
larization of the problem is, i. e., whether the result for 
the force does not change if one uses a different regular- 
ization. There are two significant points demonstrating 
such robustness. 

The first one is concerned with the geometrical regular- 
ization that we used. If one chooses domain Df. around 
qk to be not a ball but some differently shaped volume 
bounded by smooth surface Sf. ("topological ball"), then 
it is not difficult to see that all the terms of order 0(e) 
in Eq. ( p0| ) for the regularized energy remain unchanged, 
and hence our result for the force is still true. The demon- 
stration goes exactly in the same way as above, only 
the computation of the integral over the surface Sf. in 
Eq. ( |l9|) requires a well-known result from potential the- 
ory (cf. Ref. |l^). As for the first integral on the right of 
Eq. (|l7|), which defines the self-energy sei ^, its ex- 
plicit expression is not even needed, and its only relevant 
property, namely, its independence of r^, is obvious. 

An alternative way of regularization, so widely used by 
the classics during the whole "pre-Dirac delta-function" 
era, is the physical regularization, when the point charge 
qk is replaced, within the small volume D|, with some 
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smooth charge distribution of the density p e k (r) and the 
same total charge q^, and e is taken to zero in the answer. 
From the technical point of view, this approach proves 
to be more complicated in this particular case, but it 
leads again to the same O(e) terms in Eq. ( pc| ) for the 
regularized energy. The key point here is to start with 
the following expression for the regularized energy, 
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and then, instead of Eq. (g), split the potential into a sum 
of volume potentials of p k (r) over D k (which becomes sin- 
gular in the limit), and a regular addition ip e R (r). In par- 
ticular, this regularization is used by Smythe in Sec. 3.08 
of Ref. ^ for calculating the force on a single point charge 
in a domain with the zero potential at the boundary. The 
derivation there is at the 'physical level of accuracy', and 
the answer is not brought down to its physically most rel- 
evant form of Eq. (|23|). Moreover, the final answer there 
[r. h. s. of Eq. (2) in that Section] is, unfortunately, for- 
mally diverging, because of the inappropriate use of the 
notation for the total potential in place where its regular 
part should be. 

Finally, we want to end our discussion by noticing that 
the electrostatic problem we just solved, as well as its gen- 
eralizations (see Sec. |v|), involve only volume charges. On 
the other hand, magnetostatic problems that deal, for ex- 
ample, with magnetic fluxes trapped in superconducting 
media (cf. Ref. |l3|) give rise to surface charges. Analysis 
of these is of extreme importance for today's experimen- 
tal physics H]. No easy solution for the force between 
surface charges should be anticipated since the details of 
the boundary shape, such as its curvature, are expected 
to play a role; the interaction of such surface charges will 
be discussed elsewhere. 



V. GENERALIZATION: OTHER BOUNDARY 
CONDITIONS 

We can now generalize our result for other conditions 
at the boundary. A modest but potentially useful gener- 
alization is to the case of electrodes, when an arbitrary 
distribution of the potential, V(r), and not just a zero, 
is specified at the boundary: 



V(r), 
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Let us split the potential in two, 
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of which the first is caused by point charges without any 
voltage applied to the boundary, and the second is en- 
tirely due to the boundary voltage. Therefore tp^> satis- 
fies the boundary value problem of Eqs. (|l|) and (g), 
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According to what is proved above, the force on a 
charge from the field specified by the potential ip^> is 
given according to Eq. (|23|) , 
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On the other hand, potential satisfying 
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describes the field external to all the point charges, since 
it does not depend on them and their positions. There- 
fore the force exerted by this field is 



Fj 2) = -q t V^ 2 \v) 



(31) 



Using now the superposition principle, we add these two 
forces to reinstate the result of Eq. (23) in the considered 
case: 
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The mixed boundary conditions 
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where the surfaces S%, S2 are non- intersecting (Si D S2 = 
0) and comprise the whole boundary (Si U S2 = S), and 
V(r), er(r) are given functions, lead to the same standard 
result for the force [Eq. (p3|)1 without any new technical 
difficulties. Indeed, we split the total potential in two as 
in Eq. (E6J) and require that 
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The derivation of the force from ip^ goes exactly as 
in Sec. Ill and leads to Eq. (p9|). The external to the 
charges field from vp^ produces the force of Eq. d3l]), so 
by superposition the total force is again as in Eq. (|22" 
[or Eq. ©]. 
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The appropriate split of the potential in two parts 
[Eq. (p6|)1 is a little bit trickier for the Neumann bound- 
ary condition, 



dip 
dn 



= a(r), f a(r)dA + Q = 0, Q = 
Js 



N 

(38) 

Namely, the solvability criterion (the total charge must 
be zero) makes us, when splitting the potential, to add 
and subtract another charge Q (equal to the sum of the 
point charges qi) at some point r» of the domain _D, to 
obtain 



(i) = 1 
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as well as 
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(39) 
(40) 

r, r* G D- (41) 
(42) 



with both problems solvable. Again, the derivation of 
the force from ijfi-' satisfying the homogeneous boundary 
condition goes exactly as before and leads to Eq. ( |29] ) , the 
external to the charges field exerts the force given 
in Eq. (31), and by superposition the result of Eq. ( |2^ ) 
holds. The problem itself, though, is not too realistic, 
except for the case of an insulated boundary, cr(r) = 0. 
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